^ Cee 


NACA TN 3233 }ZQ 


WN ‘SAVM AHVHSITHOBL — | 


^... Sb589"UU 


I'M wet. - 


NATIONAL ADVISORY COMMITTEE 
FOR AERONAUTICS 


DUREE 


cn E 


TECHNICAL NOTE 3233 


A REVIEW OF PLANING THEORY AND EXPERIMENT WITH A 
THEORETICAL STUDY OF PURE-PLANING LIFT OF 
RECTANGULAR FLAT PLATES 
By Charles L. Shuford, Jr. 


Langley Aeronautical Laboratory 
Langley Field, Va. 


Washington 
August 1954 


p= — Lu 
P IR M ; 
TECHIES GRE T "r ^ anr 


6" uv 


LC Ny H 
WU thy fos a bt 


- —— — — -a Áo Ke 0 


TECH LIBRARY KAFB, NM 


NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS UI l | l i iil | 


| Dübb2*b 
TECHNICAL NOTE 2255 


A REVIEW OF PLANING THEORY AND EXPERIMENT WITH A 
THEORETICAL STUDY OF PURE-PLANING LIFT OF 
RECTANGULAR FLAT PLATES 


By Charles L. Sbuford, Jr. 
SUMMARY 


À summary is given of the background and present status of the pure- 
planing flat-plate lift theories. The fundamental assumptions and the 
applicability to actual calculations of the planing lift force are 
reviewed. 


A proposed theory based on the consideration of linear lifting-line 
theory less the suction component of lift plus crossflow effects is pre- 
sented. A comparison of this theory with existing planing formulas and 
experimental data is made. The agreement between the results calculated 
by the proposed theory and the experimental data is satisfactory for engi- 
neering calculations of pure-planing rectangular-flat-plate lift and cen- 
ter of pressure. 


INTRODUCTION 


Recent developments in water-based aircraft have resulted in config- 
urations utilizing planing surfaces operating in ranges of trim, length- 
beam ratio, and Froude number beyond those for which most of the avail- 
able planing theories were correlated with experimental data. In order 
to determine whether available planing theories are adequate in estimating 
the planing lift in these extended ranges, & review of these theories 
(refs. 1 to 13) and a correlation with existing data, including recent 
and unpublished data, were made and are presented herein. For purposes 
of expediency and simplification, this work is limited to the case of 
the rectangular flat plate in pure planing, that is, where buoyancy can 
be considered as negligible. 


Tn addition to this review and correlation, an additional theory 
for the 1ift and center of pressure of a rectangular flat plate was devel- 
oped and correlated with the pure-planing data. The proposed theory dis- 
tinguishes between linear and nonlinear components of lift and is divided 
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into three parts: first, a reasonably accurate approximation to the 
linear components of lift; second, an estimation of the aerodynamic 
leading~edge-suction component of lift contained in the linear term; and 
third, a method for calculating the crossflow effects. 


SYMBOLS 
A aspect ratio, b^ /s 
b beam of planing surface, ft 
Cr, lift coefficient, —E 
5 VS 
CL lift coefficient due to buoyancy, 
œ P. yes 
2 
Cr, lift coefficient based on square of beam, ze 
ib p 2 2 
— V 
2 

Cle lift coefficient based on principal wetted area, 

E Wm 

P a Lab 

2 
CN normal-force coefficient, B 

P ves 
2 

Cy speed coefficient or Froude number, V/V gb 


acceleration due to gravity, 32.2 ft/sec@ 


L lift of planing surface, 1b 
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ia 


lift due to buoyancy, lb 


wetted length of planing surface or chord of airfoil, ft 


mean wetted length, ft 


center-of-pressure location (measured forward of trailing 
edge), ft 


nondimension&l center-of-pressure location 


section lift-curve slope per radian 


normal force, lb 

free-stream dynamic pressure, Save, 1b/sq ft 

principal wetted area (bounded by trailing edge, chines, 
and heavy Spray line), sq ft 

horizontal velocity, fps 

induced vertical velocity 


induced vertical velocity at distance y from center line 
of airfoil 


distance from center line of airfoil to point where value 
of downwash is desired 


V2 Cr, 


circulation or strength of vortex, 5 


nondimensional loading parameter, Cr LS 
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A vertical load, Jb 
7 distance from center line of airfoil to vortex 
n" = mU m 
b/2 
0 = cosy 


gu = cos~ly* 
p mass density of water, slugs/cu ft 


T trim (angle between planing bottom and horizontal), radians 
unless otherwise stated 


induced angle of trim, w/V 


p: 
Ta induced angle of trim at distance y from center line of 
y airfoil, w,/V 
d 
p Pabst's aspect-ratio correction factor based on the ratio 


of wetted length to mean beam 
REVIEW OF EXISTING PLANING-LIPFT THEORY 


Wagner (ref. 1) considered the planing-force problem theoretically; 
however, his work is valuable mainly for the basic concepts presented in 
the application of the methods of airfoil theory to the plsning problem. 
Wagner's work consists of studies of the flow processes and solutions for 
the force on an idealized two-dimensional planing surface; therefore, his 
work is not directly applicable for calculating the lift on a finite- 
aspect-ratio planing surface. 


In planing theories such as that of Mayo (ref. 2) developed from 
virtual-mass considerations based on transverse flow, the assumption is 
made that the planing force can be calculated from the rate at which 
momentum is imparted to the downwash; however, the effect of aspect ratio 
(s approximated by the Pabst empirical aspect-ratio correction factor 

ref. 14). 
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In reference 35 Sokolov presented a combined theoretical and experi- 
mental solution of the planing problem. The theoretical formulas, which 
were ‘developed for the two-dimensional case, were derived by using 
Bernoulli's equation and disturbance velocities. A finite-aspect-ratio 
planing-lift formula was developed by using Sottorf's experimental results 
(ref. 15) to determine empirically the value of the factor e which is 
the ratio of the change in velocity along the planing surface to the 
velocity of the free stream. The planing formule gives the lift forces 
in three components: the hydrostatic, the one due to circulation, and 
the one due to form. The solution given by Sokolov for pure-planing 
flat-plate lift is 


CL. = e(2 - €)cos T (1) 


where curves for e are given in reference 5. 


Sokolov gives a qualitative picture of the planing problem and deter- 
mines the nature of the forces involved. The concepts presented, however, 
have not been used in the development of subsequent planing formulas, 
which have been empirical or follow the work of Wagner. 


Perring and Johnston (ref. 4) presented the empirical relationship 


Cp, = CAT (2) 


and by analyzing Sottorf's data (ref. 15) found the following formula 
to apply: 


Ops = 0.9040 2, (3) 
"Ih reference 5, Sottorf proposed the formula 

Cia = 0.8454? ‘2r ý (14) 
An equation that has a form similar to airfoil lifting-line theory 


was presented by Perelmuter (ref. 6). The equation is 


CAT 
Cr = (5) 
IS lA 
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Sedov (ref. 7) gives an equation based on the data of Sottorf 
(ref. 15) and Sambraus (ref. 16) which has the form 


_ O-7aAT 
“Ig A+ Ly (6) 


An equation that contains a linear and nonlinear term was presented 
by Siler in reference 8. The linear term was obtained by assuming a form 
similar to airfoil lifting-line theory; however, the aspect-ratio factor 
was altered to give a deflected mass one-half that predicted by Jones 
(ref. 17) for a vero-aspect-ratio wing. The nonlinear term was obtained 
by a consideration of the transverse component of the flow (see ref. 18). 
The equation can be written in the form 


C xA sin T COS T i Tc 7 


In reference 9, Korvin-Kroukovsky, Savitsky, and Lehman proposed 
an equation derived primarily on the basis of the data of Sottorf (ref. 15) 
and Sambraus (ref. 16). The formula can be written as 


Cle = 0.01249 -2(51.231) 1-1 (8) 


In reference 10, Korvin-Kroukovsky presented an equation that con- 
sisted of linear and nonlinear components. The linear term was obtained 
by a consideration of the downwash and the analytical solution for the 
potential flow about a planing surface developed by Wagner (ref. 1) and 
presented in detail by Pierson and Leshnover (ref. 19). The nonlinear 
term was obtained by a consideration of the transverse component of the 
flow. The equation has the form 


LAT 


= + 0.8812 
"Lg A+ 2 


However, this equation was empirically corrected to get better agreement 
with experimental data, so that approximately 


— 0.730 AT . 2 
Cha aa Oe (9) 
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Locke (ref. 11) assumed that the lift characteristics of low-aspect- 
ratio surfaces can be represented by a simple power function of the form 


Cla = DKr? 


where K and n depend only on aspect ratio and D is primarily a 
function of the operating conditions. For the case of the flat-plate 
planing surface Locke gives the equation 


Cis = 0.5Kra (10) 


where curves for K and n are given in reference 11. 
In reference 12 Perry assumed an equation for the ratio of planing 


lift to aerodynamic lift which converged to limits obtained by applying 
airfoil methods to the planing surface. The equation has the form 


“lg = aeons (11) 


where M represents the assumed equation for the ratio of planing lift 
to aerodynamic lift given .by 


and 


nes cos T 
l + cos T - (1 - cos diel +x sin T 
2 cos T l 


Curves for M and up are given in reference 12. 
The limit of M for zero aspect ratio is 0.88 and for infinite ` 


aspect ratio is the value p. The value of M for zero aspect ratio is 
a result given by Bollay (ref. 18). The value of M for infinite aspect 
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ratio was obtained by a consideration of the analytical solution for the 
potential flow about a planing surface developed by Wagner (ref. 1) and 
presented in detail by Pierson and Leshnover (ref. 19). 


An equation having a linear term with a form analogous to airfoil 
lifting-surface theory was proposed by P. R. Crewe of Ssunders-Roe Ltd. 
(British) in correspondence between himself and the Langley Laboratory. 
This equation, based on the data of Kaprysen and Weinstein (ref. 20), is 


Uia = gin T COS 7|2 ——L—— T 2 sin T - B sin?T (12) 
AN 2 
l + L + (3) 
where 
B= 5.0 (A > 2.0) 


Schnitzer (ref. 13) presented an equation derived from a considera- 
tion of two-dimensional deflected mass, modified for three-dimensional 
flow by the Pabst empirical aspect-ratio correction factor (ref. 14). 
The equation can be written in the form 


CL = (ER sin a cos T + 0.88 sine (13) 


PROPOSED THEORY 


An examination of experimental data indicates a pronounced nonlinear 
relationship between the planing lift coefficients and the angle of attack; 
therefore, linear theory would not provide adequate approximations to the 
planing lift. The determination of linear and nonlinear components of 
lift is the approach generally used in low-aspect-ratio airfoil theory. 
The present approach is based on the consideration of linear lifting- 
line theory less the suction component of lift plus viscous crossflow 
effects. 


AX 
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Lift 


Linear term.- The linear term is determined from & consideration of 
lifting-line airfoil theory. Since the heavy spray line (leading edge) 
of a planing surface is approximately elliptic, the airfoil theory is 
presented for an elliptic surface with elliptic loading and then modified 
for the planing case. 


By use of the Prandtl airfoil theory, the airfoil lift coefficient 
is 


Cy, = mo(T - T4) (1k) 
where mg is the slope of the section lift curve. 


If a sheet of trailing vortices located at 0.75 of the chord meas- 
ured from the trailing edge and extending to infinity behind the airfoil 
is assumed, then from vortex theory the induced angle of attack is 


b/2 
T+ = iy = a ran an (15) 
Y V AnV y/o an (n - y) 


where 1 is the distance from the center line of the airfoil to the 
vortex and y is the distance from the center line of the airfoil to 
the point where the value of downwash is desired. 


Now let 
y* = oo cos 84 (16) 
n* = m cos 0 (17) 
i = E (18) 


and 


2H 1 
pe (19) 
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y = CL1* = mgl*(v - 71) (20) 


From equation (15), 


x 
- i.f dy ao o 
Y &Jo d0 cos 0 - cos 04 


Now let 


y= D a, sin n9 (21) 


then 


2 nay, cos n0 dé 


af BS (22) 
"y Ox cos 8 - cos 0. — 01 


The solution of this equation can be obtained by means of a recurrence 
formula and the solution of & linear finite-difference equation (ref. 21); 
thus, 


V^ E in n90j (23) 
sin 04 


Since this equation is valid for any value of 04, the subscripts on T4 
and O0 can be dropped and equation (20) becomes 


Mol* 8 sin 0 


T= 2s + Y sin n9 (24) 
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For an elliptic airfoil, 


and 


Let 


then 


where 


b, = 


Let T be constant along 


bj = 


an = 


The lift coefficient 


1* = -2 sin 0 


TL 
a (T sin 0)sin nO do 
TJ 0 


the wing; then, 


T bo = bz = Bnei zc 
81 a5 = az = anl = 0 
is given by 


x. f om 
b/a 99 


M i ——À— A RG —Á— — 
- a A E mm m RS a M a eg i p, t A a Pingitore p P RI e P rg or B pa qp s PU iP HEP ap ge — 
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(25) 


(26) 


(27) 


(28) 


(29) 
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By use of equations (17), (19), and (21), equation (29) becomes 


TL 
y sin 8 dO = F Aan (30) 


-å 
Li oy 


Therefore, from equations (27) and (28), 


nc cem (31) 


which is the equation for the lift on an airfoil. 

For a low-aspect-ratio planing surface having flow only on one side, 
the lift coefficient is assumed to be one-half the value given by equa- 
tion (351) for a flat-plate airfoil, and mg, the lift-curve slope for the 


two-dimensional planing surface, is assumed to be one-half the value given 
for & flat-plate airfoil; thus, 


Cy, = O -AT (32) 


which gives the linear component of lift on a pure-planing flat plate. 
Suction component of lift.- An airfoil has a suction component of 

lift due to the large negative pressures at the leading edge of the air- 

foil; however, for a planing surface this suction component of lift does 


not appear. Therefore, the lift obtained from the linear term (eq. (22)) 
is less by an amount equal to the suction component of lift given by 


Cp œ aaae NT (33) 


which is the value indicated by Wagner in reference 22. 


The linear term (eq. (32)) less the suction component of lift is 


a O AT (ah. gine 
Oy, = S2E (n + sin?r) (34) 
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Crossflow term.- For a simple theoretical consideration of the non- 
linear term, the velocity component perpendicular to the chord is assumed 
to be of the magnitude V sin T. The drag coefficient for a planing sur- 
face of infinite aspect ratio is assumed to be 1.0, which is one-half the 
value given for a two-dimensional flat-plate airfoil. For the planing 
surface the flow is projected into components perpendicular and parallel 
to the planing-surface chord line, and the drag force associated with the 
flow perpendicular to the chord is calculated. Therefore, the normal 
force is 


= 1.0 2 S(V sin T)? 
2 
and 


CN = sin?T 


or 


è 


Cy = sin?T cos T (45) 


which is a lift due to crossflow effects and is proportional to sin?r, 
which is the concept presented for airfoils by Betz in reference 25. 


Total lift.- The total lift on & pure-planing rectangular flat plate 
can be obtained by adding equations (34) and (35) and is 


CI Ur 2r ( - sin?) + sin^T cos T (36) 
L TL A 


which represents the linear term less the suction component of lift plus 
the crossflow term. The magnitude of the crossflow effects, total lift, 
and total lift excluding suction effects is shown in figure 1. 


Comparison of proposed and previous planing formulas.- A comparison 
of the proposed theory with previous planing formules for constant length- 


beam ratios is given in figure 2. In figure 2(a) the proposed theory is 
compared with the planing formulas as presented by Sokolov (eq. (1)), 
Perring and Johnston (eq. (5)), Sottorf (eq. (4)), Perelmuter (eq. (5)), 
and Sedov (eq. (6)). In figure 2(b) the proposed theory is compared with 
the planing formulas presented by Siler (eq. (7)), Korvin-Kroukovsky 

(eq. (9)), and Schnitzer (eq. (12)). In figure 2(c) the proposed theory 


——— - - — —— — — — 
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is compared with the planing formulas presented by Korvin-Kroukovsky, 
Savitsky, and Lehman (eq. (8)), Locke (eq. (10)), and Crewe (eq. (12)). 


The values given by the formula presented by Perry (ref. 12) were 
not plotted since the results depended on the airfoil data used. Perry 
showed that by using Winter's airfoil data (ref. 24) his formula approxi- 
mated the results given by the formula presented by Korvin-Kroukovsky, 
Savitsky, and Lehman (ref. 9) for trims up to 12° and length-beam ratios 
below approximately 1.0. 


Center of Pressure 


The center of pressure on a planing surface of small aspect ratio 
may be considered to have two components, the component due to the linear- 
lift term less the suction effects and the component due to the crossflow 
term. The center of pressure for the lift due to the linear component of 
lift less the suction effects (eq. (34)) is assumed to be located at 0.75 
of the mean wetted length from the trailing edge of the planing surface. 
The center of pressure for the lift due to the crossflow term (eq. (35)) 
is assumed to be located at the center of the mean wetted length; there- 
fore, 


0.75(C - C + 0.5C 
i «OP Clota 7 Lerosstiow) * 9 Lerosstiov (37) 
Ln lade ULtotal 


which is a formula analogous to that used in airfoil theory. The com- 
ponents of lift are determined from equations (35) and (36). 


COMPARISON OF THEORY AND EXPERIMENT 


Buoyancy 


The experimental data were considered as pure planing if the lift 
coefficient due to buoyancy, calculated from the wedge-shaped volumetric 
displacement of the planing surface below the level water surface and 
given by 


. im 1 
CL. Aes "E sin 2T (38) 


NACA TN 32355 15 


did not exceed a given value. The allowable lift coefficient due to 
buoyancy, as determined from equation (38), was arbitrarily selected as 
0.01 at a trim of 16°. The allowable lift coefficient due to buoyancy for 
other trims was determined by drawing a straight line from zero trim and 
zero lift coefficient due to buoyancy through the value 0.01 at a trim 

of 16°. The permissible lift coefficient due to buoyancy for the data 
selected by this method at a trim of 2° varied from 16 percent at a length- 
beam ratio of 8 to 5.5 percent at a length-beam ratio of one-half. These 
values decreased with increasing trim so that at 30° they would vary from 
6.5 percent at a length-beam ratio of 8 to 5.1 percent at a length-beam 
ratio of one-half. Analysis of unpublished Langley tank no. 2 data showed 
that subtraction of the lift coefficient due to buoyancy computed from 
equation (38) caused the data to approximately coincide (or collapse) 

for different speeds at low trims (4°); however, at the higher trims 

(16°) only approximately one-half of this value was needed to collapse 

the data. 


Lift 


A comparison of the proposed theory with the experimental data of 
Weinstein and Kapryan (ref. 25), unpublished NACA data, data of Shoemaker 
(ref. 26), data of Locke (ref. 27), data of Sambraus (ref. 16), and data 
of Sottorf (ref. 15) is presented in figures 3 to 10. Only the experi- 
mental data indicated as pure planing by the method discussed in the 
preceding section have been considered. 


Figure 5 gives a comparison of the proposed theory with the data of 
Weinstein and Kapryan (ref. 25). Figure 4 gives a comparison of the pro- 
posed theory with unpublished Langley tank no. 2 data. The data of 
Weinstein and Kapryan were obtained for a 4-inch-beam model tested at 
various loads and speeds and the unpublished data were from a 2.5-inch- 
beam model tested at a constant speed of 20 feet per second. 


In figures 5 to 10 a comparison of experimental lift coefficients 
given in references 25, 26, 27, 16, and 15 and unpublished NACA data 
with the proposed formulas given by Crewe (eq. (12)), Locke (eq. (10)), 
Korvin-Kroukovsky, Savitsky, and Lehman (eq. (8)), and the present paper 
is presented. In general, the proposed theory gives an average of these 
data. The formula presented by Crewe (eq. (12)) is in good agreement, 
except at a trim of 30°, with the data of Weinstein and Kapryan (figs. 5(a) 
and 5(b)), with the data of Sambraus (fig. 9), and with the data of Sottorf 
(fig. 10). The formulas presented by Locke (eq. (10)) and Korvin- 
Kroukovsky, Savitsky, and Lehman (eq. (8)) do not give so good a repre- 
sentation of experimental data as the proposed theory or the formulas 
presented by Crewe (eq. (12)). There are no experimental data at high 
trims and large length-beam ratios to determine whether the planing for- 
mulas give the correct variation of lift in this region. The agreement 
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between the proposed theory and experiment is apparently satisfactory for 
engineering calculations of pure-planing rectangular-flat-—plate lift in 
the ranges where experimental data are available. 


Center of Pressure 


The variation of center-of-pressure ratio with mean wetted-length— 
beam ratio for the data of Weinstein and Kapryan (ref. 25) is shown in 
figure 11 and for the unpublished Langley tank no. 2 data in figure 12. 
The proposed theory is shown by the solid line in figures 1l and 12 where 


— = | — —. The agreement between the curve for the proposed 
b im b 

calc 
theory and the experimental points appears to be satisfactory for engi- 
neering calculations of pure-planing rectangular-flat-plate center of 
pressure in the ranges where experimental data are available. 


CONCLUDING REMARKS 


The proposed theory appears to predict the pure-planing rectangular- 
flat-plate lift and center of pressure with engineering accuracy in the 
ranges where experimental data are available; however, at high trims and 
large wetted-length—beam ratios no data are available. The correlation 
of experimental data and theory in this report seems to establish firmly 
the utility of the. singT approach (where T is trim) to the nonline- 


arity problem. 


Langley Aeronautical Laboratory, 
National Advisory Committee for Aeronautics, 
Langley Field, Va., June 10, 1954. 
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Figure 1.- Relative magnitude of components of proposed theory. 
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(a) Proposed theory and references 5 to T. 


Figure 2.- Variation of lift coefficient with angle of attack for rectangular- 
flat-plate lift formulas. 
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(b) Proposed theory and references 8, 10, and 13. 


Figure 2.- Continued. 
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(c) Proposed theory, references 9 and 11, and Crewe's equation (eq. (12)). 


Figure 2.~ Concluded. 
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Figure 5.- Comparison of proposed theory with experimental. lift coeffi- 


cients for the lh-inch-beam rectangular-flat-plate planing-surface data 
of Weinstein and Kapryan (ref. 25). 
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Figure 4.- Comparison of proposed theory with experimental lift coeffi- 
cients for a 2.5-inch-beam rectangular-flat-plate planing surface 


(unpublished NACA data). 
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Figure 7.- Comparison of proposed theory, formulas given by equations (8), 
(10), anà (12), &nà experimental lift coefficients for the 16-inch-beam 
rectangular-flat-plate planing-surface data of Shoemaker (ref. 26). 
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Figure 8.- Comparison of proposed theory, formulas given by equations (9); 
(10), and (12), and experimental lift coefficients for the 4-inch-beam 
rectangular-flat-plate planing-surface data of Locke (ref. 27). 
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Figure 9.- Comparison of proposed theory, formulas given by equations (8), 
(10), and (12), and experimental lift coefficients for the 15- and 30- 
centimeter-beam (5.91- and 11.81-inch-beam) rectangular-planing-surface 
data of Sambraus (ref. 16). 
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Figure 10.- Comparison of proposed theory, formulas given by equations (8), 
(10), and (12), and experimental lift coefficients for the 30-centimeter- 
beam (11.81-inch-beam) rectangular-planing-surface data of Sottorf 
(ref. 15). 
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Figure ll.- Variation of center-of-pressure ratio with mean wetted-Length— 
beam ratio for the l-inch-beam rectangular-flat-plate planing- surface 
data of Weinstein and Kapryan (ref. 25). 
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Figure 12.- Variation of center-of-pressure ratio with mean wetted-length— 
beam ratio for a 2.5-incb-beam rectangular-flat-plate planing surface 
(unpublished NACA data). 
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